Vector meson dominance and the 7r° transition form factor 
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It is shown that the ir° transition form factor F(Q\,Q\) differs substantially from its one-real- 
photon limit F(Qi, 0) even for rather small values of Q\ (« 0.1 GeV 2 ), which cannot be excluded in 
experiments with one "untagged" electron. It indicates that the comparison of data with theoretical 
calculations, which usually assume Q% = 0, may be untrustworthy. Our phenomenological model of 
the 7T° transition form factor is based on the vector-meson-dominance hypothesis and all its param- 
eters are fixed by using the experimental data on the decays of vector mesons. The model soundness 
is checked in the two-real-photon limit, where it provides a good parameter-free description of the 
7T° — > 2y decay rate, and in the tv° Dalitz decay. The dependence of .F(Qi,Ql) on Qi a t several 
fixed values of Q\ is presented and the comparison with existing data performed. 
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The issue of the 7r° transition form factor has recently 
attracted renewed interest in connection with the pre- 
cise measurements of the BABAR Collaboration [1], which 
seem to indicate that the asymptotic limit predicted by 
perturbative QCD [2] has been exceeded. A comprehen- 
sive review of the current theoretical situation with an 
extended list of references can be found in Ref. [3] . On 
a phenomenological side, it has recently been shown [4] 
that the vector-meson-dominance (VMD) hypothesis [5- 
7] leads to a correct description of the two-photon decay 
of the 7T° if the parameters are fixed by the data on the 
partial decay widths of the vector mesons [8]. The rate 
of this decay is related to the real-photon limit F(0,0) of 
the 7T transition form factor. It is therefore tempting to 
use the VMD also for the construction of the 7r° tran- 
sition form factor, which parametrizes the dynamics of 
the process in which two off-mass-shell photons fuse and 
form a 7r°. 

The experimental data on the tt° transition form factor 
are taken in the process e~ + e~(e + ) — > e~ +7r° + e~(e + ), 
where the photon virtualities are given by the electron 
(e + or e — ) momenta transfer squared, Qf — —q 2 , i = 1,2. 
To get the one-real-photon transition form factor only 
the data sample is utilized in which one lepton exhibits a 
small momentum transfer, e.g., \q 2 \ < 0.18 GeV 2 in the 
latest BABAR experiment [1] . 

The transition form factor F(Q 2 , Q 2 ,) is defined by its 
appearance in the 7*7*71° vertex, which describes the 
fusion of two virtual photons with the four-momenta qi 
and q2 into a tt° 



(1) 



where e is the elementary electric charge. This definition 
of the transition form factor agrees with that used in 
[1, 9], but differs from that in [10], where the factor e 2 
was absorbed in F(Qf,Q|). The two-real-photon value 
of the pion transition form factor is related to the two- 



photon decay width of the 7r° by 



r(7r°^27) = — m 3 n F 2 (0,0), 



(2) 



where a is the fine-structure constant. 

Another process in which the transition form factor 
plays a role is the decay ir° — > e + e _ 7, which was sug- 
gested by Dalitz [11] as an explanation of the anoma- 
lous events recorded in photographic emulsions exposed 
in high-flying balloons [12]. The original evaluation of 
the branching ratio 



B = 



r(7r° 



e + e 7) 



r(7r° -> 2 7 ) 



(3) 



by Dalitz as well as a later one [13], did not consider a 
possible form factor. The latter was included in [14]. The 
differential branching ratio (3) in the Berman-Geffen [14] 
variable x = M 2 /m 2 , where M is the mass of the e + e~ 
pair, reads as 



dB _ 2a (1 - x) 3 
dx Sir x 



^Wl-^ 2 (x), (4) 
x I V x 



where e = m 2 /m 2 . Form factor Fry is related to the n 
transition form factor by 



n 



F D (x) 



F(-xml,0) 
F(0,0) 



(5) 



The total branching ratio (3) is not very sensitive to 
the shape of the form factor (5) and cannot serve as a 
stringent test of theoretical calculations or phenomeno- 
logical models. Even the original Dalitz formula, which 
did not include the form factor at all, leads to B = 
1.185%, which agrees with the experimental value of 
(1.188 ± 0.035)% [15]. The form factor (5) is at small 
x usually parametrized as 



F D {x) = 1 + 



ax. 
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In 1961, Gell-Mann and Zachariasen [6] showed that the 
form factor Fjj is dominated by two resonances, namely p 
and oj, and got a positive a equal to m 2 (m p 2 +m~ 2 )/2, in 
agreement with today's observations. Until late 1980's, 
most experiments had indicated negative values of a, see 
[16] and a list of preceding experimental results there. 

Our model of the pion transition form factor F(Q 2 , Q 2 ) 
is defined in Fig. 1. In addition to the puj intermedi- 
ate state considered in [4, 6, 7], we include the following 
companions of the p°: 4>, w(1420) (denoted as oj' in what 
follows), and cj(1650) (denoted as oj"). All these reso- 
nances, except oj" ', have recently been considered in the 
model of the one-real-photon transition form factors of 
7T°, 77, and r}' [17]. 




FIG. 1. Feynman diagrams defining our model of the ir° tran- 
sition form factor. 

The combination of an isoscalar resonance with an 
isovector one in the intermediate states is unique for the 
7T° transition form factor and implies the equal contribu- 
tion of the I = and 1 = 1 currents to F(0, 0). With the 
T) and r\' , the situation is different. 

In order to evaluate the Feynman diagrams depicted 
in Fig. 1, we use the Lagrangian 



•'VpTT 



(6) 



where V is the operator of a neutral vector meson field, 
p is the isovector of the p meson fields, and tp is that 
of the pion fields. The coupling of photons to neutral 
vector mesons is given by the VMD Lagrangian (a little 
nonstandard notation of [18] adopted from [19] is used) 



= -E 



eg Vj 



in 



2 A»V U 



(7) 



where is the electromagnetic field operator. 

After comparing the amplitude corresponding to Feyn- 
man diagrams in Fig. 1 with the definition of the transi- 
tion form factor (1) we extract the latter in the following 
form 



F(QlQl) = 9 -f E G v [R v {Q\)R p {Ql) 
v=u,4>,... 

+ R v {Ql)R p {QD] , 
where Gy = Gy p7r <jy 7 and functions 

Rv{Q 2 ) 



m 2 +Q 2 



(8) 



(9) 



are the scalar parts of the vector meson propagators be- 
low the physical cut threshold in s = — Q 2 , sth — ~m 2 . 



Tensor parts do not contribute thanks to the presence of 
the Levi-Civita tensor in (6). For R p (Q 2 ) we will alter- 
natively use the form 



R P (Q 2 ) 



M 2 (0) 



m 2 (-q 2 ) + q 2 



(10) 



with the running mass squared given by the dispersion 
formula 

00 

M 2 (s) = M 2 (0) - S -Vf ffi&dS (11) 

p 1 TT J S'(S' — S) 



and satisfying M 2 (m 2 



i 2 and dM 2 /ds = at s = 



to 2 . The energy dependent total width T p includes the 
contributions from the following final states: 7r + 7r~, KK, 
7t°cj, ijnn, 7r°7, 777, and ir + ir~-f. For details, see [20]. 

Inserting (8) into (5) and calculating the derivative of 
Fd(x) at x = we are getting the following expression 
for the slope of the Dalitz decay form factor 



rat 



a = 



M 2 (0) 



1 \ - Gy 

J2 v G v ^m 2 v 



(12) 



Here, D is the derivative of the running mass (11) at 
s = 0. When a fixed mass form (9) is used also for 
R P (Q 2 ), then D = and M 2 (0) = to 2 . If we, in addition, 
keep only the first term V = lo in the sums over V above, 
we recover the Gell-Mann-Zachariasen [6] formula. 

The values of coupling constants Gy p7r are determined 
as follows. For the vector mesons with masses below or 
close to the pir threshold [w(782), 0(1020)] we consider 
the decay chain V — > p + n — > 3ir. The comparison 
of the decay width formula given in [18] with the rec- 
ommended value from [15] yields the product Gy p7r g pl 
where g p is the coupling constant of the usual pirn inter- 
action Lagrangian. Using the value g 2 = 35.70 ± 0.19, 
as it follows from the p meson decay width [15], we 
get Gy p ^. For the vector mesons with higher masses 
[a/, oj" , J/ip(lS)] it is sufficient to explore the simpler 
V — > p + 7r decay width formula. Concerning the oj', 
the Review of Particle Physics [15] gives only intervals 
for its mass and total width and no quantitative esti- 
mate for the pir branching fraction. We use therefore 
the values [21] to^ = (1.38 ± 0.02 ± 0.07) GcV/c 2 and 
IV = (0.13±0.05±0.10) GeV as measured by the BaBAR 
Collaboration [22] and B(oj' -> p + it) = (69.9 ± 2.9)% 
from the wavelet analysis [23] of the e + e~ annihilation 
data. In the case of the oj", we again use the BABAR [22] 
values m^, = (1.667 ± 0.013 ± 0.006) GeV/c 2 and T^,, = 
(0.222 ±0.025 ±0.020) GeV with the wavelet analysis [23] 
branching fraction B(w" -> p + tt) = (38.0 ± 1.4)%. The 
resulting values of the coupling constants Gy plx squared 
are summarized in Table I. The coupling constants G U(OT 
and G^p-x themselves differ in sign, as it follows from the 
analysis of the p — > decay [18] and from the SU(3) 
symmetry [19]. This results in the negative sign of G^ 
shown in Table I. 
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TABLE I. The squares of the coupling constants in La- 
grangians (6) and (7) obtained from the vector meson decay 
data described in the text. Also shown are the parameters 
Gv, which enter the transition form factor (8). 



V 


Gv P Tv (GeV- 2 ) 


Ov-y X 10 2 


G v (GeV" 1 ) 


w(782) 


216.2 ±3.0 


1.375 ±0.046 


1.724 ± 0.031 


0(1020) 


0.676 ± 0.020 


2.214 ± 0.031 


-0.122 ±0.002 


w(1420) 


11.7 ±1.1 


0.20 ±0.17 


0.152 ±0.136 


w(1650) 


3.97 ±0.61 


0.76 ±0.11 


0.174 ±0.026 



Now, we determine the coupling constants in the VMD 
Lagrangian (7). The value g% = 4/g 2 p = 0.1120(10) fol- 
lows from the normalization of the charged pion form fac- 
tor. The squares of other coupling constants gv-y are eval- 
uated from the dilepton decay width of the corresponding 
vector mesons. In the case of the lo, the e + e~ decay width 
from [15] is used. For the (/>, we obtain it as a product of 



the full width and the e+e~ 



<7„» 7 , we combine the full width from 



branching fraction. To get 
22] (shown above) 



with the i 



branching fraction B(lo'' 



) = (32± 



1) x 10 7 from [23]. Let us note that the product of the 
pir and e + e~ branching fractions from [23] agrees with 
B(uj" -> e+e-)B(uj" -> p + w) = (1.3 ±0.1 ±0.1) x 10~ 6 
got later by the BaBAR Collaboration [24]. The same 
concerns the full width. Unfortunately, the situation 
with the lo' is more controversial. The product of the 
pir and e + e~ branching fractions from the wavelet anal- 
ysis [23] is about twofold of that observed by BABAR [24] . 
In order not to overestimate the contribution of the lo' 
to the 7T° transition form factor, we use the BABAR value 
B(uj' -> e + e-)B(co' -> p+w) = (0.82±0.05±0.06) x 10~ 6 
to get G 2 , = G^^g 2 ^, from which the estimate of g 2 ,^ 
presented in Table I is obtained. 

As all the parameters of our model of the ir° transi- 
tion form factor are now determined, we are ready to 
test its soundness by evaluating the characteristics of the 
two-photon and Dalitz decays of the neutral pion and 
comparing them to the measured values. To avoid the 
cumulation of the errors, we do not use Table I for the 
error analysis, but calculate the errors of the results di- 
rectly from the input data (masses, full widths, branch- 
ing fractions, and their errors). To account for possi- 
ble correlations among the input quantities, we sum the 
contributions to the final errors from various sources lin- 
early. Two of the calculated quantities, namely the mean 
lifetime of ir° and the Dalitz decay slope parameter are 
shown in Table II for various versions of our model, i.e., 
for various choices of the isoscalar vector mesons enter- 
ing the sum in Eq. (8) and for two possible treatments 
of the p propagator. The branching ratio of the Dalitz 
decay to the two-photon decay is not shown in the table. 
It acquires the same value of 1.196(1)% in all versions 
of the model, in agreement with the experimental value 
of 1.188(35)% [15]. Our results concerning the 7r° tran- 
sition form factor are presented in Figs. 2 and 3. The 



TABLE II. Comparison of various versions of our model with 
the 7T° decay data. When converting the calculated ty" — > 
2-/ decay rate into the tv° mean lifetime r, the experimental 
branching fraction of 0.988 was used. The results above the 
horizontal divider were obtained with the constant p mass, 
those below it with the running mass (11). 



V 


r x 10 17 (s) 


a x 10 2 


LO 
L0,(j> 


7.6 ±0.4 
8.8 ±0.6 


3.002 ± 0.002 
3.049 ± 0.004 


U,tj> 
U),(j>,Lo' 
LO, (f>, Oj" 
LO,(j>,Lo',Lo" 


8.8 ±0.6 
7.3 ±1.9 
7.2 ±0.7 
6.1 ± 1.7 


3.282 ± 0.003 
3.190 ±0.097 
3.163 ±0.024 
3.089 ± 0.099 


Data [15] 


8.4 ±0.5 


3.2 ±0.4 



contributions of the various n°p°V vertexes to the one- 
real-, one- virtual-photon form factor together with their 
sum are shown in Fig. 2. They were all calculated us- 
ing the running mass (11) of the p meson. In addition, 
the sum of the same contributions but calculated using 
the fixed p mass is depicted by a dotted curve. It differs 
from the running-mass case only marginally. The strong 
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FIG. 2. The individual contributions to the pion transition 
form factor _F(Q 2 ,0) multiplied by Q 2 and their sum (solid 
curve) with the uncertainties originating in the errors of input 
parameters. The sum of individual contributions (not shown) 
that were calculated assuming the constant mass of the p(770) 
is depicted by the dotted curve. 

dependence of the transition form factor F(Q 2 ,Q 2 ) on 
the virtuality Q\, which is given by the four-momentum 
transfer of the "untagged" electron, is demonstrated in 
Fig. 3. The expression Q\F{Q\,Q^) is presented as a 
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FIG. 3. Dependence of QlF(Ql,Ql) on Q{ for four fixed 
values of Q\ calculated in the full version of our model. The 
CELLO [10], CLEO [9], and BABAR [1] data are also shown. 



function of Q\ for four different virtualities Q\ from 
to 0.3 GeV 2 . The full version (V = oj,<f>,w', and lu") of 
our model with running mass of the p is used. Compari- 
son with existing data indicates that our model would be 
able to describe them if the (unmeasured, but certainly 
nonvanishing) absolute value of the momentum transfer 
squared of the "untagged" electron decreased with the 
rising Q\. A more quantitative account is given in Ta- 
ble III, which shows that the agreement of our model with 
the data below Q\ — 9 GeV 2 is excellent but deteriorat- 
ing for higher Q\. There are two possible ways of improv- 
ing our model in an attempt to get a better description of 
the high-Q 2 data. First, the inclusion of higher isovector 



[p(1450), ,9(1700)] and isoscalar [0(1680), 0(2170)] res- 
onances. Second, a consequent use of the running-mass 
propagators also for other resonances [not only for p(770) 
as here] , which guarantee the correct analytic properties 
of the transition form factor. 



TABLE III. The comparison of our parameter-free model with 
the data in various Q\ ranges. The virtuality Q 2 of the photon 
radiated from the "untagged" electron was assumed the same 
for all n data points within a particular Q% range and deter- 
mined by minimizing the x 2 ■ The confidence levels (C.L.) are 
also shown. 



Qi (GeV 2 ) 


n 


x 2 


Ql (GeV 2 ) 


C.L. (%) 


0-9 


27 


10.67 


0.187 ±0.008 


99.7 


9-18 


6 


3.77 


0.107 ±0.013 


58.3 


18-36 


4 


4.47 


0.010 ±0.025 


21.5 



The effect of nonvanishing Q\ has already been quanti- 
tatively studied in terms of the photon momentum asym- 
metry parameter A = (Q\ — Qi)/(Ql + Qi) within the 
spectral quark model [25] and the incomplete vector- 
meson dominance (IVMD) [3]. In the former, the transi- 
tion form factor decreases with rising Q 2 (falling A)-see 
Fig. 2 in [25]-as in our model. In the latter, the tendency 
is inverse-see Fig. 7 in [3]. This difference is caused by 
not keeping the IVMD parameters c and My constant, 
but allowing them to vary in order to get the best fit for 
each particular A. 

Our model, in spite of its deficiencies, supports the 
conclusion of Refs. [3, 25] that it is important to pay 
more attention to the dependence of the transition form 
factor on both virtualities in theoretical calculations and 
experimental analyses. 
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